Casimir-Polder interaction between an atom and a dielectric grating 
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We develop the scattering approach to calculate the exact dispersive Casimir-Polder potential 
between a ground-state atom and a rectangular grating. Our formalism allows, in principle, for 
arbitrary values of the grating amplitude and period, and of the atom-grating distance. We compute 
numerically the potential for a Rb atom on top of a Si grating and compare the results with the 
potential for a flat surface taken at the local atom-surface distance (proximity force approximation). 
Except for very short separation distances, the potential is nearly sinusoidal along the direction 
transverse to the grooves. 



I. INTRODUCTION 

Vacuum field fluctuations are modified close to mate- 
rial surfaces, resulting in the usually attractive Casimir- 
Polder force pQ on a nearby ground-state atom. This 
effect has been measured by a number of experimen- 
tal techniques, including deflection of atomic beams [2J, 
classical [31 13] and quantum reflection by the attractive 
Casimir-Polder potential [SHZ], and dipole oscillations of 
a Bose-Einstein condensate (BEC) close to a dielectric 
surface [5]. 

More recently, experiments involving non-trivial ge- 
ometries paved the way to new applications. The mea- 
sured reflection probability of a BEC from a Si surface 
with a square array of closely spaced thin pillars increased 
by a factor of nearly four as compared to a flat sur- 
face [9]. A Si surface with an array of wall-like parallel 
ridges was shown to work as a reflection diffraction grat- 
ing for atoms incident along a nearly grazing direction 
[lOj . Several diffraction orders in the quantum reflection 
from a microstructured grating consisting of Cr strips 
on a flat quartz substrate were measured [llj. Alterna- 
tively, nano-fabricated transmission atom gratings allow 
for a direct measurement of the dispersive atom-surface 
potential in the non-retarded van der Waals regime [12] , 

Remarkable experimental progress has also been 
achieved in the closely related field of Casimir interac- 
tions between material surfaces |13j . Recent experiments 
have revealed interesting geometry effects in the normal 
force between a Si rectangular (lamellar) grating and a 
metallic spherical surface [14j [15] and in the lateral force 
between two metallic gratings [16j . 

This ensemble of new experiments clearly motivates 
the theoretical analysis of how geometry molds the quan- 
tum field fluctuations giving rise to the Casimir-Polder 
and Casimir interactions. In this paper, we develop a 
non-perturbative theory for the Casimir-Polder interac- 
tion between a ground state atom and a dielectric grating. 
Our results are valid, in principle, for arbitrary values of 
the grating amplitude. Preliminary results were already 



applied in the analysis of quantum vortex generation in a 
BEC induced by the Casimir-Polder potential of a rotat- 
ing grating |17) . Here we present a detailed derivation of 
the Casimir-Poder potential and a variety of numerical 
examples that illustrate its main characteristics. 

Calculations beyond the simple planar geometry are 
extremely involved because different frequency and 
length scales contribute to the interaction. Since disper- 
sive forces are not additive, it is not possible to build 
up the atom-surface potential from the more elemen- 
tary atom-atom interaction. In the pairwise summation 
(PWS) approach, the non-additivity is corrected by a 
'calibration' provided by the planar case at a given sepa- 
ration distance range |18j . Since the non-additivity cor- 
rection is geometry and distance dependent, PWS results 
are barely more accurate [19 than the results obtained by 
taking the proximity force approximation (PFA) |20j . in 
which the potential is simply approximated by the planar 
case taken at the local atom-surface distance. 

Only recently new theoretical tools have been devel- 
oped for analyzing non-trivial geometries |21j beyond the 
PFA and PWS approaches. In the scattering approach 
[2"2l [25] . the interaction potential is written in terms 
of reflection operators describing non-specular scattering 
(diffraction) by non-planar surfaces. This allows for a 
description of the Casimir effect that captures its full ge- 
ometry dependence. In the specific case of the Casimir- 
Polder interaction between a ground-state atom and a 
material surface, the exact potential can be calculated for 
arbitrary distances, taking due account of finite response 
times of both atom and material medium, provided the 
surface reflection operator is known |24j . 

By computing the reflection operator of the grating as 
a perturbation of the planar symmetry, Refs. [2H [25] de- 
rived the Casimir-Polder potential to first order in the 
grating amplitude. The resulting expressions are valid 
only when the grating amplitude is the smallest length 
scale in the problem. However, in order to enhance the 
non-trivial geometry effects associated with the depar- 
ture from the planar symmetry, it is of course interesting 
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to have large amplitudes. This was indeed the case in 
the atom-surface [10] and surface- surface [14] [16] ex- 
periments carried out so far. In this paper, we compute 
the exact reflection operator by applying the differen- 
tial theory of diffraction gratings , and then deriving 
the exact Casimir-Polder potential from the scattering 
approach [24] . Our approach is similar to the theory 
developed in Refs. [2"7rf2T)] for treating surface-surface in- 
teractions. 

Non-perturbative results were previously obtained for 
a toy model describing the Casimir interaction of a scalar 
field satisfying ideal Dirichlct boundary conditions with 
a small sphere above a corrugated surface [3U]- In this 
paper, we develop a full electromagnetic theory that takes 
into account the electromagnetic responses of real atoms 
and material surfaces (we consider a Rb atom above a 
Si grating in the numerical examples). This allows us to 
cover the entire range of separation distances, from the 
unretarded short-distance van der Waals regime to the 
Casimir-Polder large-distance asymptotic limit. 

The paper is organized as follows. In Sec. II, we derive 
the formal results for the Casimir energy which are ap- 



plied in the numerical implementations discussed in Sec. 
III. Sec. IV presents our concluding remarks. 



II. SCATTERING APPROACH TO THE 
ATOM-GRATING INTERACTION 

We consider a spherically symmetric ground-state 
atom located at (xa,Ua,za) above a non-planar surface 
(see Fig. [T]), corresponding to a profile function h(x,y) 
giving the local surface height with respect to a refer- 
ence plane at z = 0. It is convenient to develop the 
scattering formula in the plane- wave basis |k, ±,p) [22j . 
where k is the two-dimensional wave-vector component 
parallel to the xy plane, +(— ) represents upward (down- 
ward) propagation direction and p stands for polariza- 
tion. We assume that the atom-surface separation dis- 
tance za — h(xA,y A ) is much larger than the atomic di- 
mensions, allowing us to approximate the potential to 
first order in the atomic polarizability a(uj) or, equiva- 
lently, to first order in the atomic reflection operator [21] 



(k,-,p\K A mV,+,p') = -f a ^e i ( k - k ')-^e-("+"')^e;(k,iO ^(k'.O, (1) 

with ta — (xa,Ua) an d k = \fk 2 + £ 2 /c 2 representing the wave-vector z-component associated with the imaginary 
frequency £ (k' is defined in terms of k' in the same way), (k, i£) are unit vectors corresponding to a given 
polarization basis (to be chosen later in this section). 




FIG. 1: Atom on top of a rectangular grating. 

The zero-temperature Casimir energy is then obtained by expanding the general scattering formula [22 , 23 to first 
order in TZa- This corresponds to neglecting all multiple reflections between the atom and the surface (reflection 
operator TZs) but the single round-trip containing one reflection by each one (see Ref. [24] for a detailed derivation). 
The Casimir potential is then written as an integral over the positive imaginary semi-axis in the complex frequency 
plane: 

r°° H£ f H 2 k f H 2 k' 

u(x A , yA ,z A ) = -hl g I I ^J2(^+Mnsm\k\-,p')^\-,p'\n A (iO\K+,P)- (2) 
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An exact expression for the Casimir-Polder interaction face can also be written in terms of an integral over real 
energy between a ground-state atom and a generic sur- 
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frequencies. Once such an integral is computed with an 
appropriate choice of integration contour in the complex 
frequency plane, one recovers Eq. ^ , which is in terms of 
purely imaginary frequencies uj = i£. On one hand, the 
real frequency formalism suffers from highly oscillatory 
contributions, but allows for the identification of the rel- 
evant modes (evanescent and propagating) of the electro- 
magnetic field that contribute to the interaction energy. 
For example, both modes are known to play a fundamen- 
tal role in the Casimir energy in the plane-plane geometry 
[3"T] . We expect that both kinds of modes also have a key 
role in the Casimir-Polder atom-grating interaction (this 
analysis is beyond the scope of the present paper). On 
the other hand, the formalism in terms of purely imag- 
inary frequencies has a smooth, exponentially decaying 
frequency behavior, but it does not allow for the identifi- 
cation of the separate contributions from evanescent and 
propagating modes. However, Eq. ^ does include the 
total contribution of all modes. 

In this paper, we calculate the exact surface reflec- 
tion operator 1Z$ for the rectangular (lamellar) grating 
shown in Fig. [I] The periodic profile function h(x) has 
period d, amplitude a and corresponds to a groove width 
s. The reference plane z = is located at the bottom of 
the groove. Translational symmetry along the y direc- 
tion implies that 1Z$ does not change k y . When defining 
the polarization basis, we exploit this symmetry as in 
waveguide theory: polarization H is such that E y = 
with E x , E z , H x and H z given in terms of H y by using 
Maxwell equations. For the electric field components at 
a real frequency w propagating in a medium of dielectric 
constant e, we have E x — —iajjd z H y /(euj 2 — c 2 k y ) and 
Ef = icud x H y /(euj 2 — c 2 k 2 ). In vacuum (e = 1), the 
corresponding unit vectors are then given by 



for some integer j and k x ^ G [— ir/d, ir/d]. Diffraction by 
the periodic grating will give rise to new Fourier compo- 
nents modulated by integer multiples of 2ir/d. Although 
these new Fourier components correspond to different in- 
tegers f (with f — j representing a given diffraction or- 
der), they all correspond to the the same k x ^ in the first 
Brillouin zone according to Bloch's theorem. For a given 

incident plane wave \(k x ^ + j 2 f)it + k y y, —,p), we add, 
in the homogenous region above the grating (z > a), a 

superposition of reflected plane waves |(fci°^ + f 2 f)x.+ 
k y y,+,p') (Rayleigh expansion), with amplitudes given 
by the matrix elements (f \p'\1Zs\j,p) of the grating re- 
flection operator. In the bulk region z < 0, the dielec- 
tric constant is also uniform and the field is written as 
a simple plane wave expansion in terms of transmission 
amplitudes. On the other hand, within the inhomoge- 
neous grating region < z < a, we derive the non-trivial 
z dependence of the field by solving coupled differential 
equations (differential approach) [53]. We then solve for 
the matrix elements (j' ,p'\lZs\j,p} by matching the dif- 
ferent field expansions at z = a and z = 0. 

The properties of the reflection operator 71$ discussed 
above allow us to simplify the general expression ^ for 
the Casimir-Polder potential energy. By plugging the 
result ([lj for the atomic reflection operator into ([2]), we 
find 



h r dg 
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Polarization E is likewise defined by the condition H y 



0, with J5| = 
ic 2 k y d z H y /(euj 2 
vectors are given by 



ic rkyd x Ey/ (eto 



c 2 k 2 ) and Ef 



kz), so that the corresponding unit 



Juj 2 /c 2 
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^— k x x =F \/lu 2 /c 2 — k 2 : 



-^Jlo 2 /c 2 



(4) 



Since the surface profile is periodic along the x direc- 
tion, the fields are pseudo-periodic functions of x (Bloch's 
theorem): E(x + d, y, z) = e lfe ^° >d E(x, y, z) for some k x ^ 
in the first Brillouin zone [— n/d, Tt/d]. The wave- vector 
x-component of a given incident plane wave can always 
be cast in the form 



k^ = k^+j 



2n 



(5) 



The sums over the Brillouin zones j and j' run from — oo 
to oo, the polarizations p,p' are either E or H and Kj — 

\J ki 3 '' 2 + k 2 + £ 2 /c 2 [see Eq. (5)]. The scalar products 
are calculated from the expression for the unit vectors 
given by Eqs. ([3| and Q. When replacing 
yjufljt?- k 



it, 



we find 



= 1 



[.t')j.„.„ ., 

e/c 2 + k 2 y ' 

r 2 k 2 



C r^ y \r^x ry>j ">x rvj 

ae+c 2 ki) 



The potential ^ does not depend on y A and is a 
periodic function of xa (with the same period of the 
grating) as expected. By choosing the origin x = 
at the mid-point of one of the grooves, the matrix el- 
ements {j' ,p'\R-s\jiP} turn out to be real and consis- 
tent with the even parity of the surface profile. They 
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must also satisfy reciprocity relations [35] which read, 
for our polarization basis, Kj(j,p\1Zs(kx , k y , ,p') = 
(26 pp < - l)K j >{-j',p'\ns(-k ( °\-kyj£)\- j,p). By com- 
bining these properties, we may cast ^ as U(xa, za) — 
Y IJ ,j> C jJ'( z A)e 27TiU ~ ] ' )xA/d in terms of real coeffi- 
cients satisfying Cjjt(zA) — Cj',j(zA)- Thus, Eq. ^ 
yields a real potential satisfying the required relations 
U(xa,za) — U(—xa,za) = U(d — xa,za) (reflection 
symmetries with respect to the groove and plateau mid- 
points). 



III. NUMERICAL RESULTS 

In this section, we present numerical results for the 
Casimir-Polder potential given by (|6|. We consider a 
Rb atom, whose dynamic polarizability a(z£) is obtained 
from Ref. [33] . interacting with an intrinsic silicon grat- 
ing. The corresponding dielectric constant e(i£), required 
for the evaluation of the grating reflection operator IZs, 
is computed from available data at real frequencies [34] 
by using a Kramers-Kronig relation [35] . 

When evaluating ^ , we need to truncate the sum over 
the Brillouin zones j, j' at some finite value iV max , so that 
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the number of zones is 2N„ 



1. The value of A max re- 



quired for a given accuracy decreases with increasing dis- 
tance za, because the factor e~ tyKi+K i'^ ZA kills the contri- 
bution of large values of j and j'. By comparing results 
obtained from different values of iV max , we found that 
]V max = 3 was sufficient to achieve an accuracy at the 
level of a percent for the numerical examples considered 
below. 

It is instructive to compare our exact results with 
those obtained within the proximity force approximation 
(PFA) [35] , which here corresponds to computing the po- 
tential U for a given surface profile from the potential 
for a planar surface taken at the local atom-surface 
distance: U(x a ,Va,z a ) ~ U^[z A - h(x A ,VA)]- For 
the rectangular grating considered here, deviations from 
PFA are quantified by the ratio p = U(xa, za) /U^ (za) 
when the atom is on top of a groove and by p — 
U(xa,za)/U^(za — a) when the atom is on top of a 
plateau. 

In Fig. [2] we plot p as a function of distance for an atom 
above the plateau midpoint xa = d/2 and several differ- 
ent combinations of amplitude and period, with groove 
width s = d/2. We always find p < 1, in qualitative agree- 
ment with PWS, with the grooves making the potential 
less attractive than the potential of a homogeneous plane 
surface at za = a. When the atom is very far from the 
surface, za a, it feels the corrugation as a very small 
perturbation of the plane symmetry and then the poten- 
tial approaches the value for a planar surface (p = 1). 
As the distance decreases, the departure from the plane 
geometry becomes increasingly important. 

The curves corresponding to the same ratio a/d 
are close but do not coincide, as expected since non- 
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FIG. 2: Variation of p = U(xa, Za)/U(°' {za — a) versus za/o 
and (za — a)/d (inset) for an atom above the plateau midpoint 
xa = d/2. We take fixed values of a and d for each set of data 
points. The angular aperture 9 = 2 arctan[d/4(z j 4 — a)] shown in 
the diagram controls the departure from the PFA limit. 



geometric length scales associated to characteristic fre- 
quencies of the Rb atom and the Si bulk are also relevant. 
For large distances, the dominant field frequencies are 
much smaller than the atom and medium characteristic 
frequencies, and then the potential is dominated by the 
instantaneous response associated to the zero-frequency 
polarizability and dielectric constant. Thus, the potential 
depends only on the ratios between the geometric lengths 
a, d and za in this asymptotic limit. On the other hand, 
as the distance decreases, finite response times of both 
atom and medium give rise to a richer scenario. The 
crossover between the long-distance Casimir-Polder and 
the unretarded van der Waals regimes is in the range 
200 — 300 nm for a Rb atom interacting with a Si surface 
(see for instance Fig. 4 of Ref. [24] ). When the separa- 
tion distance decreases below za/cl = 5, it already cor- 
responds to the transition to the van der Waals regime 
for a = 50 nm but not yet for a = 100 nm. The values for 
a = 50 nm then move away from the curve for a = 100 nm 
and the same a/d = 1/6 as shown in Fig. p] A similar 
effect, for a = 100 nm, 200 nm and a/d = 1/3 is also ap- 
parent in Fig. [2] 

As p approaches the PFA limit p —> 1 at shorter dis- 
tances, the amplitude a is no longer the most relevant 
length scale capturing the variation with distance. In 
the inset of Fig. [2j we plot p as a function of (za — a)/d, 
which is directly related to the angular aperture 9 = 
2 arctan[d/4(zA— a)] of the plateau width as seen from the 
atom location above the plateau midpoint (see diagram 
in Fig. [2}. When 9 is very close to 180° [{z A -a)/d < 1], 
we expect the effect of the border of the plateau to be 
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negligible, and then the potential should be well approx- 
imated by the result for an infinite plane (PFA). The 
inset of Fig. [2] shows that the variable (za — a)/d indeed 
captures the main effect behind the departure from PFA 
as the atom is displaced away from the surface, since the 
curves corresponding to different values of amplitude and 
period collapse near each other at short distances. For 
all different values of a and d shown in the figure, the 
point of maximum 1 — p is near (za — a)/d ~ 0.8, corre- 
sponding to an angular aperture 9 ss 35°. The maximum 
1 — p depends mainly on a/d and is larger for larger a/d 
as expected. 
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FIG. 3: Variation of p = U(xa, za)/U(°) (za) versus z^/a for an 
atom above the groove midpoint xa = 0. 

When the atom is above the groove region, the poten- 
tial is stronger than the potential for a planar surface 
at the bottom of the groove (p > I), again in quali- 
tative agreement with the PWS picture. In Fig. [3j we 
plot p as a function of distance for x = (groove mid- 
point) and s = d/2. The curves corresponding to differ- 
ent (fixed) values of a and d merge at distances za/cl ~ 8. 
This remarkable property was first found in Ref. [30] in 
the context of the scalar Casimir-Polder model. It may 
be interpreted as an effect of averaging the small length 
scales associated to the surface profile when considering 
the long-wavelength field fluctuations that provide the 
main contribution at large distances. 

At shorter distances, the curves move apart following 
the order of larger a/d. Although atomic/material fre- 
quency scales play a relevant role at za < 300 nm, the 
values for different a and d with the same a/d shown in 
Fig. [3] are still relatively close. 

From an overall comparison of Figs. [2] and [3j we may 
conclude that the relative deviation from PFA \p — 1| is 
generally larger above the groove than above the plateau. 
In Fig. |4j we plot \p — 1| as a function of the lateral 
position for a fixed local distance za — H(xa) = 3a, with 




0.1 0.2 0.3 0.4 0.5 



\ z_ 

d/2 x 

FIG. 4: Lateral variation of \p — 1| at a fixed local distance za — 
h(xA) with a = 100 nm and d = 600 nm. We take za = 3a above 
the groove (0 < xa < d/2), and za = 4a above the plateau {d/2 < 
xa < d). 

a = 100 nm, d = 600 nm, and s — 300 nm. PFA is indeed 
worse above the groove, where the deviation increases 
rapidly as one approaches the corner. Above the plateau, 
PFA is also worse near the corner, but the variation of p 
with lateral position is much smaller. 

At a given plane z — za, the potential profile along 
the x axis is approximately sinusoidal for the relatively 
large distances considered here. In Fig. [5^, we compare 
the values for U(xa, za)/U^(za — o,/2) with a sinusoidal 
variation (solid line) corresponding to period d and am- 
plitude [U{d/2,z A ) - U(0,z A )]/2. We use z A = 3a and 
the same parameters employed in Fig. [4] 

At very large distances, the main contribution comes 
from specular reflection at the first Brillouin zone j = 
j' — 0, yielding a flat potential according to Eq. ^ 
(see Figs. [2]and[3|. As the distance decreases, the ma- 
jor correction, corresponding to the first diffraction or- 
der j — j' — ±1, produces a sinusoidal modulation with 
period d. Deviation from the sinusoidal shape (coming 
from \j — j'\ > 2) is such that the lateral force stiffness 
\d 2 U/dx A \ is slightly larger at the groove midpoint than 
at the plateau midpoint, as shown in Fig. [5^,. 

The deviation from the sinusoidal shape and the am- 
plitude of oscillation are controlled by the angular aper- 
ture parameter (za — a)/d, whereas ZA/a controls the 
convergence of the spatial average (U(xa, za)) to its 
large distance limit J7 (0) (z A - a/2) = [C/ (0) (^ - a) + 
U tS)) (z A )]/2 + 0(a/z A ) 2 when s = d/2. For instance, in 
Fig. [5|d we plot the potential as a function of XA/d for 
Za = d = 2a = 200 nm. These parameters correspond to 
Za/cl — 2 and (za — a)/d = 0.5, to be compared with 
z A/ a — 3 and (za — a)/d — 0.33 for Fig. [5^,. While 
the average potential is closer to the asymptotic value 
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FIG. 5: Lateral variation of the Casimir-Polder potential divided 
by U(°\za ~ a/2) for a = 100 nm with (a) d = 2z A = 6a; and (b) 
d = za = 2a. The solid lines represent the sinusoidal function of 
period d that coincides with the exact potential at the groove and 
plateau midpoints. 

U^\za — a/2) in Fig. [H^, (larger za/o), the amplitude 
is smaller and the shape closer to the sinusoidal fit in 
Fig. [5]d (larger (za — a)/d). As the atom is displaced 
away from the surface, the potential becomes more si- 
nusoidal and with a smaller amplitude. With the same 
parameters of Fig. [SJd, the amplitude decreases by two 
and three orders of magnitude at za — 3a and za = 4a, 
respectively. 



IV. CONCLUDING REMARKS 

We have shown that PFA overestimates the potential 
strength above the plateau. As the atom is displaced 
away from surface, the deviation 1 — p is mainly con- 
trolled by the parameter (za — a)/d associated to the an- 
gular aperture of the plateau as seen by the atom. The 
maximum deviation takes place near (za — a)/d = 0.8 
and is larger for deeper grooves (larger a/d) as expected. 
As the distance is increased past this value, the potential 
again approaches the result for a planar surface. On the 
other hand, the potential is underestimated by the PFA 
above the groove. The deviation increases strongly as the 
atom is laterally displaced closer to a corner, where the 
presence of the ridge walls have a stronger effect. 

The naive PWS picture of the Casimir-Polder inter- 
action agrees qualitatively with several features found 
in this paper. Nevertheless, PWS typically underes- 
timates the magnitude of the Casimir-Polder potential 
in the grating geometry [19] . particularly in the non- 
perturbative regime considered here |17| . For instance, it 
predicts an anisotropy of the Casimir-Polder potential far 
too small to nucleate vortices in a Bose-Einstein conden- 
sate by rotation of the grating, whereas the exact theory 
presented here predicts an anisotropy above the required 
values under realistic experimental conditions [17] . 

We have shown that the contribution of higher diffrac- 
tion orders to the potential is increasingly small as the 
distance increases, resulting in a sinusoidal variation 
along the transverse x direction for (za — a)/d ~ 0.5 
(aperture angle smaller than 53°). 

The scattering approach developed in this paper 
might be readily adapted to consider more general two- 
dimensional periodic patterns, like the pillars structure 
employed in the quantum reflection experiment reported 
by Ref. [5]. Recently, the Casimir interaction between 
two material surfaces imprinted with two-dimensional pe- 
riodic patterns has been analyzed in detail |37j . The 
method developed here thus paves the way for the quanti- 
tative analysis of a variety of interesting geometries with 
applications in enhanced atomic quantum reflection and 
diffraction. 
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